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Search results for nucleon decays p e + X, p —>• p + X, n —¥ v 7 (where X is an invisible, massless 
particle) as well as dinucleon decays np —> e + v, np —>• p + v and np —> t + v in the Super-Kamiokande 
experiment are presented. Using single-ring data from an exposure of 273.4 kton ■ years, a search 
for these decays yields a result consistent with no signal. Accordingly, lower limits on the partial 
lifetimes of r p ^ e + x > 7.9 x 10 32 years, r p ^ p + x > 4.1 x 10 32 years, T n .^ vl > 5.5 x 10 32 years, 
r np _ >e + I/ > 2.6 x 10 32 years, r np _^ li + v > 2.2 x 10 32 years and r 7lp _ >r + J/ > 2.9 x 10 31 years at a 90% 
confidence level are obtained. Some of these searches are novel. 

PACS numbers: 12.10.Dm,13.30.-a,11.30.Fs,14.20.Dh,29.40.Ka 


Many signs indicate that the Standard Model (SM) 
of particle physics is an incomplete description of na¬ 
ture. Gauge coupling unification, charge quantization 
and other features suggest a more unified account, such 
as a Grand Unified Theory (GUT) 0-S , as an under¬ 
lying fundamental theory. While the unification scale 
(~ 10 15 — 10 16 GeV) is unreachable by accelerators, rare 
processes predicted by these theories, such as proton 
decay, can be probed by large underground detectors. 
Being a signature prediction of GUTs, observation of 
an unstable proton would constitute robust evidence for 
physics beyond the SM, while non-observation will strin¬ 
gently constrain theoretical models. 


The simplest unification scenarios based on minimal 
SU(5) and supersymmetric (SUSY) SU(5) have been de¬ 
cisively ruled out by bounds on p -} e + 7T° 0-0 and 
p —> vK + |. Many alternative scenarios as well as po¬ 
tential signatures are possible (see 0 ] for review). 

In this Letter, we analyze a broad class of nucleon 
and dinucleon decay channels with a showering or non¬ 
showering single Cherenkov ring signature within the 
Super-Kamiokande (SK) experiment, using the technique 
of spectral fit First, we have considered two gen¬ 

eral two-body decays p —> e + X and p —>- p + X , where X 
is a single unknown invisible particle which is assumed to 
be massless. These searches are distinct from the model- 
dependent inclusive analyses of H, I3| listed in the PDG 
U|. Similarly, we also consider n —> ^ 7 . Though this 


radiative process is suppressed, it has a clean signature 
and has been considered in the context of SU(5) [15|, with 
some models [ 9 ] predicting a lifetime of 10 38±1 years. 

While single nucleon A B = 1 processes have been in 
general well studied, dinucleon A B = 2 channels also 
pose great interest. These higher-dimensional processes 
can become significant in models which suppressproton 
decay and could be connected to baryogenesis [l6|, ac¬ 
counting for the observed baryon asymmetry of the uni- 


171 ]. Such a connection may already be hinted at 


verse 

from the requirement of baryon number violation as a 
necessary condition for explaining the asymmetry [l 8 j . 
The disappearance A B = 2 reactions, with invisible fi¬ 
nal state particles, have been studied and no signal ex¬ 
cess was observed lj|-21|. The channels np —> e + v, 
np — » p + is and np — » t + v violate baryon number by 
two units and violate lepton number by either two or 
zero units. They can become significant in models with 


an extended Higgs sector [lfl Izl. which could be con¬ 
sidered in the context of GUTs [23( . While the r cannot 
occur in single nucleon decay, in dinucleon decay the r 
channel is allowed [24^ . The process np —» t + v has not 
been experimentally studied before and in addition to the 
electron and muon channel searches we present the first 
search in the r channel. 


In this work, SK data is analyzed from an exposure of 
the 22.5 kton fiducial mass for 273.4 kton-years, covering 
four running periods (SK-I through SK-IV). Details of 
the detector design and performance in each SK period, 


can be found in |25|, [26|]. This analysis considers only 
events in which all observed Cherenkov light was fully 
contained (FC) within the inner detector. 

Since final-state neutrinos or X (by definition) are not 
observed, the only signature of p —> e + X, p —» p + X, 
n —>• 1 / 7 , np —>■ e + v and np —> p + u is a single charged 
e + or pX lepton, or single 7 . Thus, the invariant mass 
of the initial state cannot be reconstructed and the sig¬ 
nal will be superimposed on a substantial atmospheric 
neutrino background in the e-like and /i-like momentum 
spectra. For the np —> t + v decay, only the r —> e + vv 
and r —>• p + vv channels are considered, with the respec¬ 
tive branching ratios of 17.8% and 17.4%. This allows 
us to perform all the analyses within a unified frame¬ 
work. The previous searches for n —> wy, np —>• e + v and 
np —> p + v, which were performed with a smaller detector 
using a counting method, resulted in the lifetime limits 
of 2.8 x 10 31 years 0 , 2.8 x 10 3 ° years H 3 and 1.6 x 10 3 ° 
years [10, respectively. In contrast, the spectral fit em¬ 
ployed within this work, allows utilization of the extra 
information from the energy dependence of signal, back¬ 
ground and the systematic errors. This methodology has 
been recently employed in the SK nucleon decay Uj, 27 1 
and dark matter analyses j28| . 

The nucleon decay signal events are obtained from 
Monte Carlo (MC) simulations, in which all the nucleons 
of the water H 2 O molecule are assumed to decay with 
equal probability. The final state particles are gener¬ 
ated with energy and momentum uniformly distributed 
in phase space. The effects of Fermi motion, nuclear 
binding energy as well as nucleon - nucleon correlated de¬ 
cays [ 29 ] are taken into account for both nucleon H3-S2 
and dinucleon searches [331 ]. The signal Fermi momentum 
distributions are simulated using a spectral function fit to 
electron- 12 C scattering data [34|. The SK detector sim- 
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ulation [26[ is based on the GEANT-3 [35| package, with 
the TAUOLA [36] package employed for decaying the r 
leptons. For the np —> t + v mode we generated three 
MC samples, with the r decaying to e + vv : p + vv and all 
decay channels. The latter allows us to study sample con¬ 
tamination in the two selected leptonic r channels from 
the hadronic r channels and thus identify sample purity 
after the event selection. We have confirmed that the re¬ 
sulting MC charged lepton spectra from r —> e + vv and 
r —>■ p + vv decays agree with the theoretical formula [37ij. 
Forp —>• e + X andp —>• p + X modes, the invisible X parti¬ 
cle cannot be a fermion by spin conservation, but in our 
spin-insensitive MC it was simulated as a neutrino. In 
total, around 4,200 signal events were generated within 
the fiducial volume (FV) for each SK period for single 
nucleon decays and around 8,400 for dinucleon decays. 

Atmospheric neutrino background interactions were 
generated using the flux of Honda et. al. [38| and the 
NEUT simulation package [39j, which uses a relativistic 
Fermi gas model. Background MC corresponding to a 
500-year exposure of the detector was simulated for each 
detector phase. We used the same atmospheric neutrino 
MC as the standard SK oscillation analysis [40]. 

The event selection applied to the fully-contained data 
is the following: (A) a single Cherenkov ring is present, 
(B) the ring is showering (electron-like) for p —> e + X, 
n — > V"/, np —>• e + v and np —> t + v (t —>• e + vv) and 
non-showering (muon-like) for p —> /x + 2f, np p + v and 
np — > t + v (t — > p + vv), (C) there are zero decay elec¬ 
trons for modes with an e-likc ring and one decay elec¬ 
tron for those with a pi-like ring, (D) the reconstructed 
momentum lies in the range 100 MeV/c < p e < 1000 
MeV/c for p —>• e + X , n — > v'y and in the range 200 
MeV/c < p M < 1000 MeV/c for p —»• p + V, with the 
range extended to 100-1500 MeV/c for dinucleon decays 
with an e-like ring and 200-1500 MeV/c for those with 
a p-like ring. In total, approximately 37,000 FC events 
were obtained in the SK-I to SK-IV data-taking periods. 
After the criteria (A)-(D) have been applied, the final 
data samples for single nucleon decay searches with an 
e-like ring contain 8,500 events and 6,000 events for the 
case of pi-like ring, with momenta up to one GeV/c. To 
search for dinucleon decays we consider lepton momenta 
up to 1500 MeV/c. The final samples for the dinucleon 
modes contain 9,500 events for the e-like channels and 


6,500 events for the pi-like ones. See Ref. [4l| for details 
regarding reconstruction. 

The signal detection efficiency is defined as the frac¬ 
tion of events passing selection criteria compared to the 
total number of events generated within the true fidu¬ 
cial volume. The average detection efficiency for e-like 
channels is 94.0 ± 0.4% for all SK data-taking periods. 
For the /x-like channels, the average detection efficiency 
is 76.4 ± 0.6% for SK-I to SK-III and 91.7 ± 0.4% for 
SK-IV. The increase in efficiency observed in SK-IV for 
channels with a /x-like ring, comes from a 20 % improve¬ 


ment in the detection of muon decay electrons after an 
upgrade of the detector electronics [261 ]. 

For the e-like momentum spectrum up to 1500 MeV/c, 
the dominant background contribution, composing 75.8% 
of the events, comes from the v e charged-current (CC) 
quasi-elastic (QE) neutrino channel. The v e CC single- 
pion production constitutes 13.0% of the background, 
while the z/ e CC coherent-pion, CC multi-pion and 
neutral-current (NC) single-pion productions contribute 
around 1 . 1 %, 1 . 1 % and 1.6%, respectively. About 3.5% 
and 1.1% of events come from zy, NC single-pion and 
coherent-pion production. For the pt-like momentum 
spectrum up to 1500 MeV/c, the dominant contribution 
of around 78.6% comes from CCQE. Similarly, CC 
single-pion, CC coherent-pion and CC multi-pion as well 
as NC single-pion production contribute around 16.2%, 
1.4%, 1.6% and 0.8%, respectively. 

After event selection, a spectral fit is performed on 
the reconstructed charged lepton momentum distribution 
of the events. The y 2 minimization fit is based on the 
Poisson distribution, with the systematic uncertainties 
accounted for by quadratic penalties (“pull terms”) 0 - 
The y 2 function used in the analysis is 


nbins 

y 2 = 2 ^ ( N^ + N\ 
2=1 

AT exp = 


obs 


ATobs 

ln-Isp -1 
N ex p 


iVs, 


V C-f 

' rr ■ 


.7=1 


XV a , 


jyback _|_ p . 7Vj Sig | (1 + ^ f’f), 

1=1 


(1) 


where i labels the analysis bin. The terms N° hs , N^ 1S , 
V l back , V® xp are the numbers of observed data, signal 
MC, background MC 

and the total (signal and background) MC events in 
each bin i. The index j labels the systematic errors, while 
ej and fl correspond to the fit error parameter and the 
fractional change in the 7V® xp bin due to 1-sigma error 
uncertainty (Jj , respectively. The fit is performed for two 
parameters a and fj, which denote the background and 
signal normalizations, respectively. After the event se¬ 
lection, the signal MC distribution is normalized to the 
background by the integral, which in turn is normalized 
to the SK livetime. This allows us to identify the fit point 
(a,/3) = (1,0) with the no-signal hypothesis. Similarly, 
(a, /3) = (0,1) signifies that the data is described by sig¬ 
nal only, with the signal amount equal to background MC 
normalized (pre-fit) to livetime. The y 2 minimization is 
carried out over each a and /3 in the grid according to 
d\ 2 /dej = 0. The resulting global minimum is defined 
as the best fit. Further details on the fit and specifics of 
systematic error treatment can be found in [Til, 28, 42|. 
For the np —> t + vv mode, after the appropriate event se¬ 
lection is applied to both MC samples of r —> e + vv and 
r —> /x+i/zq the samples are combined for the fit, allowing 
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TABLE I: Systematic errors of spectrum fits, with la uncertainties and resulting fit pull terms. Errors specific to 
signal and background are denoted by S and B, while those that are common to both by SB. 


Decay mode 


p —> e + X 

p —> p + X 


Systematic error 

1-a uncertainty (%) 

Fit pull {a) 

Fit pull (a) 


Final state interactions (FSI) 

10 

0.10 

-0.60 

B 

Flux normalization (E v < 1 GeV) 

25 a 

-0.23 

-0.08 

B 

Flux normalization (E v > 1 GeV) 

15 b 

-1.44 

-0.50 

B 

Ma in v interactions 

10 

0.69 

0.23 

B 

Single meson cross-section in v interactions 

10 

-0.55 

-0.14 

B 

Energy calibration of SK-I, -II, -III, -IV 

1.1, 1.7, 2.7, 2.3 

0.58, -0.91, 0.48, 0.38 

-0.54, 0.07, -0.14, 0.26 

SB 

Fermi model comparison 

10 c 

-0.08 

0.70 

S 

Nucleon-nucleon correlated decay 

100 

0.00 

0.06 

S 


a Uncertainty linearly decreases with logE^ from 25% (0.1 GeV) to 7% (1 GeV). 

b Uncertainty is 7% up to 10 GeV, linearly increases with log E„ from 7% (10 GeV) to 12% (100 GeV) and then 20% (1 TeV). 
c Estimated from comparison of spectral function and Fermi gas model. 


TABLE II: Best fit (a,/3) parameter values, best fit y 2 / d.o.f., no signal Ay 2 , 90% C.L. value of /3 parameter, 
allowed number of nucleon decay events in the full 273.4 kton ■ years exposure and a partial lifetime limit for each 
decay mode at 90% C.L. The sensitivity and lifetime limit for dinucleon decay modes are per 16 0 nucleus. 


Decay mode 

Best fit 

Best fit 

No signal 

Data 

Data 

Sensitivity 

t/B 


(ol,P) 

y 2 /d.o.f. 

Ay 2 

/3gocL 

A^ocl 

(xl 0 3i yr.) 

(xl 0 3i yr.) 

p —>■ e + X 

(1.050, 0.002) 

70.9/70 

0.19 

0.013 

108 

79 

79 

n —> 1/7 

(1.045, 0.004) 

70.5/70 

0.43 

0.015 

125 

58 

55 

p —> p + X 

(0.960, 0.016) 

63.2/62 

3.43 

0.032 

187 

77 

41 

np —> e + v 

(0.955, 0.000) 

122.5/110 

0.00 

0.004 

33 

10 

26 

np —> p, + v 

(0.910, 0.000) 

97.0/102 

0.00 

0.005 

36 

11 

20 

np —>■ r + i/ 

(0.910, 0.000) 

224.6/214 

0.00 

0.006 

96 

1 

3 


us to obtain a single value for the permitted number of 
nucleon decays at 90% CL. 

The systematic errors can be divided into signal- 
specific (S), background-specific (B) as well as detector 
and reconstruction errors, which are common to both 
signal and background (SB). The two signal specific sys- 
tematics are from Fermi motion and nucleon - nucleon 
correlated decay. For background, in order to method¬ 
ically select the dominant systematics, we started from 
more than 150 errors employed in the SK oscillation anal¬ 
ysis 0 and chose those which affect the analyses bins by 
more than 5% (|//| > 0.05). Relaxing this criteria to 1% 
does not significantly alter the results, but complicates 
the analysis [11} . As in [11} , we have found that the dom¬ 
inant contributions originate from uncertainties related 
to neutrino flux and energy calibration (common to both 
signal and background). Including the signal systemat¬ 
ics, the total number of considered errors is 11 and they 
are the same for all modes. In Table Q] we display the 
complete list of systematics, their uncertainties and fit¬ 
ted pull terms for two representative examples p —> e + X 
and p —>• p + X. 

The spectral fit determines the overall background and 
signal normalizations a and /3, with the fit results dis¬ 
played in Table QI] The outcome shows that no signif¬ 


icant signal excess has been observed, with the data 
Ay 2 = y 2 — Xrnin being within lcr of the background 
only hypothesis for all search modes except for p —> p + X, 
which is within 2 cr. 

The lower lifetime limit on the processes can then 
be computed from the 90% confidence level value of /3 
(/^ 90 Cl)? which translates into the allowed amount of 
signal at 90% confidence level according to iVgocL = 
/^ 9 ocl ' A^ SIgnal , where V slsnal is the total number of signal 
events. The partial lifetime limit is then calculated from 


rgoCh/13 


E SK4 \ /\rnucleons 

sk—SKI A sk ' e sk • W _ 

A90CL 


( 2 ) 


where B is the branching ratio of a process, e s k and A s k 
are the signal efficiency and the exposure in kton • years 
for each SK phase, A 90 cl is the amount of signal allowed 
at the 90% confidence level and j\r nucleons is the number of 
nucleons per kiloton of water, corresponding to 3.3 x 10 32 , 
2.7x 10 32 and 3.3x 10 31 for proton, neutron and dinucleon 
decay searches, respectively. 

The resulting fitted spectra for the 273.4 kton • years 
of combined SK data can be found in Figure [1] The 
upper figures display best-fit result for atmospheric neu¬ 
trino background (solid line) without signal fitted to data 
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FIG. 1: [Top] Reconstructed momentum distribution for 273.4 kton • years of combined SK data (black dots) and 
the best-fit result for the atmospheric neutrino background Monte Carlo (solid line). The corresponding residuals 
are shown below, after fitted background subtraction from data. [Bottom] The 90% confidence level allowed nucleon 
decay signal (hatched histograms), from the signal and background MC fit to data. All modes are shown (overlaid), 
with e-like channels on the left and p-likc channels on the right. 


(black dots) and the corresponding residuals after the fit¬ 
ted MC is subtracted from data. It is seen that the back¬ 
ground MC describes the data well. The bottom figures 
display the 90% C.L. allowed signal (hatched histogram), 
obtained from the fit of background with signal to data, 
with all the e-like and p,-like spectra overlaid with all the 
modes. The Nqqcl as well as resulting sensitivities and 
calculated lifetime limits for the decays are shown in Ta¬ 
ble EO The sensitivities were obtained assuming that data 
are described by background. For the np — >■ tv mode we 
have combined the r channels e + vv and p + vv, weighted 
by their respective branching ratios. This limit is then 
multiplied by 1.15 to account for roughly 85% sample 
purity of the tau channels. We set the lower limits on 
the partial lifetimes of the decay modes at the 90% C.L., 
with the results shown in Table HT1 

In conclusion, the single Cherenkov ring momentum 
spectra in Super-Kamiokande are well described by at¬ 
mospheric neutrinos, including the effect of neutrino os¬ 
cillation and systematic uncertainties, up to 1500 MeV/c. 
We find no evidence for any contribution from the six 
different nucleon and dinucleon decay modes that would 


produce a showering or non-showering Cherenkov ring. 
The results of this analysis provide a stringent test of 
new physics. The obtained limits represent more than 
an order of magnitude improvement over the previous 
analyses of n —> vj [Bj and two orders of magnitude 
for np —> e + v and np —> p + v The searches for 

p —> e + X, p —> p, + X (where X is an invisible, massless 
particle) and np —> t + v are novel. The dinucleon de¬ 
cay limits restrict A B = 2 processes with L violated by 
either zero or two units. 
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tion, Culture, Sports, Science and Technology, the U.S. 
Department of Energy, and the U.S. National Science 
Foundation. 


* Deceased. 

[1] J. C. Pati and A. Salam, Phys. Rev. D8, 1240 (1973). 



























6 


[2] J. C. Pati and A. Salam, Phys. Rev. Lett. 31, 661 (1973). 

[3] H. Georgi and S. Glashow, Phys.Rev.Lett. 32, 438 
(1974). 

[4] H. Fritzsch and P. Minkowski, Annals Phys. 93, 193 
(1975). 

[5] C. McGrew, R. Becker-Szendy, C. Bratton, J. Breault, 
D. Cady, et al. (IMB-3 Collaboration), Phys.Rev. D59, 
052004 (1999). 

[6] K. Hirata et al. (KAMIOKANDE-II Collaboration), 
Phys.Lett. B220, 308 (1989). 

[7] M. Shiozawa et al. (Super-Kamiokande Collaboration), 
Phys.Rev.Lett. 81, 3319 (1998), hep-ex/9806014. 

[8] K. Kobayashi et al. (Super-Kamiokande Collaboration), 
Phys.Rev. D72, 052007 (2005), hep-ex/0502026. 

[9] P. Nath and P. Fileviez Perez, Phys.Rept. 441, 191 
(2007), hep-ph/0601023. 

[10] G. Fogli, E. Lisi, A. Marrone, D. Montanino, and 
A. Palazzo, Phys.Rev. D66, 053010 (2002), hep- 
ph/0206162. 

[11] V. Takhistov et al. (Super-Kamiokande Collaboration), 
Phys.Rev.Lett. 113, 101801 (2014), 1409.1947. 

[12] J. Learned, F. Reines, and A. Soni, Phys.Rev.Lett. 43, 
907 (1979). 

[13] M. Cherry, M. Deakyne, K. Lande, C. Lee, R. Steinberg, 
et al., Phys.Rev.Lett. 47, 1507 (1981). 

[14] K. Olive et al. (Particle Data Group), Chin.Phys. C38, 
090001 (2014). 

[15] D. Silverman and A. Soni, Phys.Lett. B100, 131 (1981). 

[16] J. M. Arnold, B. Fornal, and M. B. Wise, Phys.Rev. D88, 
035009 (2013), 1304.6119. 

[17] L. Canetti, M. Drewes, and M. Shaposhnikov, New 
J.Phys. 14, 095012 (2012), 1204.4186. 

[18] A. Sakharov, Pisma Zh.Eksp.Teor.Fiz. 5, 32 (1967). 

[19] C. Berger et al. (Frejus Collaboration), Phys.Lett. B269, 
227 (1991). 

[20] T. Araki et al. (KamLAND Collaboration), 
Phys.Rev.Lett. 96, 101802 (2006), hep-ex/0512059. 

[21] M. Litos, K. Abe, Y. Hayato, T. Iida, M. Ikeda, et al., 
Phys.Rev.Lett. 112, 131803 (2014). 

[22] R. N. Mohapatra and G. Senjanovic, Phys.Rev.Lett. 49, 


7 (1982). 

[23] L. Arnellos and W. J. Marciano, Phys.Rev.Lett. 48, 1708 
(1982). 

[24] D. Bryman, Phys.Lett. B733, 190 (2014), 1404.7776. 

[25] Y. Fukuda et al. (Super-Kamiokande Collaboration), 
Nucl.Instrum.Meth. A501, 418 (2003). 

[26] K. Abe et al., Nucl.Instrum.Meth. A737, 253 (2014), 
1307.0162. 

[27] K. Abe et al. (Super-Kamiokande Collaboration) (2013), 
1305.4391. 

[28] K. Choi et al. (Super-Kamiokande), Phys.Rev.Lett. 114, 
141301 (2015), 1503.04858. 

[29] T. Yamazaki and Y. Akaishi, Phys.Lett. B453, 1 (2000). 

[30] H. Nishino et al. (Super-Kamiokande), Phys.Rev. D85, 
112001 (2012), 1203.4030. 

[31] C. Regis et al. (Super-Kamiokande Collaboration), 
Phys.Rev. D86, 012006 (2012), 1205.6538. 

[32] K. Abe et al. (Super-Kamiokande Collaboration), 
Phys.Rev. D90, 072005 (2014), 1408.1195. 

[33] J. Gustafson et al. (Super-Kamiokande), Phys.Rev. D91, 
072009 (2015), 1504.01041. 

[34] K. Nakamura, S. Hiramatsu, T. Kamae, H. Muramatsu, 
N. Izutsu, et al., Nucl.Phys. A268, 381 (1976). 

[35] R. Brun, F. Carminati, and S. Giani (1994). 

[36] S. Jadach, Z. Was, R. Decker, and J. H. Kuhn, Corn- 
put. Phys.Commun. 76, 361 (1993). 

[37] M.-C. Chen and V. Takhistov, Phys.Rev. D89, 095003 
(2014), 1402.7360. 

[38] M. Honda, T. Kajita, K. Kasahara, S. Midorikawa, 
and T. Sanuki, Phys.Rev. D75. 043006 (2007), astro- 
ph/0611418. 

[39] Y. Hayato, Nucl.Phys.Proc.Suppl. 112, 171 (2002). 

[40] K. Abe et al. (Super-Kamiokande), Phys. Rev. D91, 
052019 (2015), 1410.2008. 

[41] M. Shiozawa (Super-Kamiokande Collaboration), 
Nucl.Instrum.Meth. A433, 240 (1999). 

[42] R. Wendell et al. (Super-Kamiokande Collaboration), 
Phys.Rev. D81, 092004 (2010), 1002.3471. 


